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Abstract. We present the spectral analysis of an XMM-Newton observation of the X-ray binary 4U 1700+24 , 
performed during an outburst in August 2002. The EPIC-PN spectrum above 1 keV can be modeled by a blackbody 
plus Comptonization model, as in previous observations. At lower energies, however, we detect a prominent soft 
excess, which we model with a broad Gaussian centered at ~ 0.5 keV. In the high resolution RGS spectrum we 
detect a single emission line, centered at 19.19^q'q| A. We discuss two possible interpretations for this line: O VIII 
at redshift z = 0.012±g;Sg| or Ne IX at redshift z ~ 0.4. 
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1. Introduction 

Based on its optical identification with the M2 III giant 
■ HD154791 (Garcia et al. 1983; Masetti et al. 2002; also 
known as V934 Her), the X-ray source 4U 1700+24 is 
classified as both a Low Mass X-ray Binary (Liu et al. 
2001) and a symbiotic- like binary (Garcia et al. 1983). The 
lack of signatures of binarity and the initially marginal 
positional coincidence (Morgan & Garcia 2001) originally 
made this identification uncertain. More recently, the de- 
tection of radial velocity variations from the proposed op- 
tical counterpart with the same period as X-ray brightness 
variations (Galloway et al. 2002) has added support to the 
identification. 

In X-rays, 4U 1700+24 has long periods in which it is 
faint, and occasional episodes in which the flux increases 
by more than an order of magnitude. Assuming a dis- 
tance of 420 pc (Masetti et al. 2002), the 2-10 keV lumi- 
nosity of 4U 1700+24 varies from 2xl0 32 , which is much 
higher than expected from an isolated M giant of this spec- 
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tral type, to 10 34 erg s" 1 . Masetti et al. (2002) reported 
the analysis of six observations from five different X-ray 
satellites, including an observation with the Rossi X-Ray 
Timing Explorer (RXTE; Jahoda et al. 1996) taken dur- 
ing a ~ 100-day outburst in November 1997. Although no 
coherent pulsations or QPOs were detected, their analysis 
shows that the X-ray flux is also variable on short time 
scales (~ seconds). 

In July 2002, the RXTE All-Sky Monitor light curve 
of 4U 1700+24 indicated that the source was undergoing 
a new outburst, and a Target of Opportunity observation 
was performed by the XMM-Newton satellite. In the fol- 
lowing sections we discuss the analysis of the spectroscopic 
X-ray data from this observation. 

2. Observations and data reduction 

4U 1700+24 was observed with the XMM-Newton satellite 
(Jansen et al. 2001) on 2002 August 11. The XMM-Newton 
X-ray optics consist of three nested Wolter-I mirror as- 
semblies, illuminating five X-ray detectors which operate 
simultaneously. CCDs sensitive to photons in the range 
0.15-15 keV are situated at the focal point of each the 
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three mirrors, two MOS and one PN types, comprising 
the European Photon Imaging Camera (EPIC; Striider et 
al. 2001, Turner et al. 2001). Two of the mirrors also il- 
luminate the Reflection Grating Spectrometers (RGS; den 
Herder et al. 2001) which disperse photons in the range 
5-35 A onto a pair of dedicated off-axis CCDs. The PN 
and MOS cameras have spectral resolutions of ~80 eV and 
~70 eV (at 1 keV) and point spread functions (PSFs) of 
~6" and ~5" for the PN and MOS cameras, respectively. 
The RGS spectral resolution is 0.04 A although it has a 
significantly smaller effective area (~100 cm 2 at 1 keV for 
each of the two units, compared to ~1200 cm 2 for the PN 
and ~900 cm 2 for the sum of the two MOS detectors) and 
cannot produce 2-dimensional images. 

The exposure time of the 4U 1700+24 observation 
(corrected for the dead-time) was 6 ks for the PN, 7 ks for 
the MOS, and 8 ks for the RGS. To avoid excessive optical 
loading on the PN and MOS cameras from the bright op- 
tical companion of 4U 1700+24 (V=7.6), the "thick" filter 
was selected. The X-ray brightness of 4U 1700+24 during 
this high-state observation necessitated the use of Small 
Window operating mode, in which only a small portion 
(about 4'x4' for the PN and 2'x2' for the MOS) of the 
CCD is read out in a reduced frame time (6 ms for PN 
and 300 ms for MOS). The short readout time of this mode 
helps reduce errors due to photon pile-up, which occurs 
when two (or more) photons hit nearby pixels during a 
single CCD exposure, producing an event that is indistin- 
guishable from a single, higher energy photon. The main 
effects of photon pile-up are an incorrect reconstruction 
of the event energies and an underestimation of the source 
count rate. 

All the data were processed using the XMM-Newton 
Science Analysis System (SAS, version 6.1.0) 1 . The use of 
Small Window mode only partially mitigated the effects 
of pile-up 2 in the PN and MOS data, and so in our ex- 
traction of source photons we excluded the central part of 
the source PSF, where the incident count rate per pixel is 
the highest. The EPIC spectra were thus extracted from 
annular regions centered on the source, with inner and 
outer radii of 10" and 40", respectively. With these re- 
gions, we were able to completely remove the portion of 
the image affected by pile-up from the PN data, but not 
from the MOS data. The choice of an even larger inner 
radius for MOS would have greatly reduced the number 
of valid counts. Therefore, the MOS spectra are not dis- 
cussed in this paper. To generate the PN spectrum, only 
single and double events were retained, and the resulting 
spectrum was rebinned both to have at least 30 counts 
in every energy channel and at most three bins per PN 
energy resolution element. A background spectrum was 
extracted from a region far enough from 4U 1700+24 to 
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2 The presence of pile-up in a dataset can be investigated by 
examining the energy distribution of single, double, and mul- 
tiple events (composed, respectively, of one, two, and greater 
than two adjacent pixels being above the detection threshold) 



avoid contamination by source photons. In the 0.2-12 keV 
range, the background contributes less than 0.5% of the 
total count rate for most of the observation, increasing to 
3% during the last 1500 seconds. 

The RGS spectra (source and background) were ex- 
tracted, and response matrices calculated using the stan- 
dard reduction procedures. Since the source photons are 
widely spread along the dispersion axis, the RGS data are 
not affected by pile-up. The RGS source spectra were also 
rebinned in order to have at least 30 counts in every chan- 
nel. 

3. Analysis and results 

3.1. Spectral Continuum 

In order to compare the XMM-Newton spectrum of 
4U 1700+24 with previous observations, we fitted the PN 
spectrum in the 0.3-12 keV range with the two-component 
model used by Masetti et al. (2002): a blackbody plus 
Comptonization (COMPST; Sunyaev & Titarchuk 1980), 
with absorption fixed to the Galactic value of N H = 
4xl0 20 cm~ 2 (Dickey & Lockman 1990). The resulting 
fit is unacceptable (x 2 ed = 4.05 for 273 degrees of free- 
dom) due to a large excess at low energies. The fit can be 
improved by adding a component that contributes below 
~ 1 keV. We obtain the best results by modeling this soft 
excess with a broad (a ~0.1 keV) Gaussian line centered 
at ~0.5 keV. We report the spectral parameters in Table 
1. We show the PN spectrum and its residuals with re- 
spect to this model in Figure 1. Although even our final 
fit is only marginally acceptable (Xred = 1-19 for 270 d.o.f., 
corresponding to a null hypothesis probability of 0.018), 
the addition of a systematic error of 1.5%, which is com- 
patible with the current calibration uncertainties of the 
PN, reduces xl e d to 1. No narrow features are detected in 
the PN spectrum. In particular, the 3er upper limit on the 
equivalent width of a Fe K-a fluorescent line at 6.4 keV is 
~ 15 eV. 

3.2. The emission line at 19.2 A 

The high spectral resolution of the RGS allows us to 
search for narrow features in the X-ray spectrum of 
4U 1700+24 at low energies. In Figure 2, we show the 
first-order RGS spectrum of 4U 1700+24 and residuals 
with respect to the best-fit model derived from analysis 
of the PN data. In addition to some deviations from the 
continuum model due to problems in the cross-calibration 
between the PN and RGS instruments (Kirsch et al. 2004) , 
an emission line is apparent at about 19 A (Fig. 2). The 
addition of a Gaussian component in the model accounts 
for the residuals in this region of the spectrum (see bottom 
^pa^icl of figure 2) and improves the x 2 from 1.65 (for 765 
cfo.f.) to 1.54 (for 762 d.o.f.). The emission-line parame- 
ters are: A = 19.19±g;g| A, a = 3.9±f J cV, line flux equals 
(4.9+i'g) x 10~ 4 photons cm~ 2 s _1 , and equivalent width 
equals 9.7+ 3 j eV (all the quoted errors are 3cr errors). 
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Table 1. Best fit parameters of the PN spectrum in the 0.3-12 keV energy range. Errors are at the 3 a confidence 
level for a single interesting parameter. 



X// (10- < m 
kT e 



*zrr 

(keV) 
T 



kT BB (keV) 
Rbb (m) (ft) 

En nc (eV) 

a (eV) 
Equivalent width (eV) 



4 (fixed) 

18.5±J;5 
0.0229±0.0007 
1.40±0.03 
113±3 
505^ 



(a) Normalization of the COMPST model in XSPEC version 11.3 
^ Radius at infinity assuming a distance of 420 pc 
Absorbed flux in units of ergs cm~ 2 s _1 



- + 13 

-12 



1 1 9+ 25 



198 



+34 
28 



Xted (d.O.f.) 


1.188 (270) 


0.3-2 keV flux (c) 


8.4xl0~ n 


2-10 keV flux (c) 


5.3xl0~ 1() 




Energy (keV) 

Fig. 1. The PN count rate spectrum of 4U 1700+24 , and 
residuals (in units of a) with respect to the best-fit model 
described in Table 1. The individual model components 
are also shown: dashed line for the broad Gaussian, dot- 
ted line for the blackbody and dash-dotted line for the 
Comptonization model. 



In the X-ray atomic lines catalogue ATOMDB 1.3.1 3 , 
only very low emissivity transitions are found within 3<r 
of the line position; most are satellite lines of H-like and 
He-like oxygen, with a maximum emissivity of 8.4xl0 -18 
photons cm s _1 . The Ly-e line of N VII (which is not 
included in the ATOMDB database) has a wavelength 
compatible with the observed emission line (A = 19.118 
A, Verner et al. 1996a), but the probability for this tran- 
sition is very low. The Ly-cv. transition of H-like oxygen 
(O VIII) has an emissivity more than 2 orders of magni- 
tude higher than that of the lines mentioned above (which 
makes it typically one of the strongest lines found in cos- 
mic X-ray sources) and a rest-frame wavelength of 18.97 
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A, slightly smaller than the value we find for the emission 
line shown in Figure 2. It is therefore possible that the 
X-ray emission line in 4U 1700+24 is O VIII at redshift 
z — u.uiz_ 004 . 

The accuracy of the RGS wavelength scale is better 
than 10 mA (Pollock 2004), with an additional 2.3 mA of 
systematic error for each arcsecond error in the pointing or 
source coordinates. For the 4U 1700+24 observation, the 
additional error is smaller than 10 mA, since the EPIC 
positional accuracy is better than 4" and the RGS data 
were processed using source coordinates derived from the 
EPIC images. Therefore, we exclude the possibility that 
the wavelength shift from the O VIII rest-frame position 
is due to an incorrect wavelength scale in the RGS instru- 
ment. 

Other smaller structures are also present in the resid- 
uals. In particular, the data differ substantially from the 
model around the instrumental oxygen edge at ~ 23.5 
A, where the RGS effective area calibration is complex. 
We can improve the RGS spectral fit slightly ( ~x^ ed =1.5l 
for 761 d.o.f.) by including an overabundance of neutral 
oxygen in the photoelectric absorption model (with cross 
section from Verner at al. 1996b). Keeping the hydrogen 
column density fixed at the value of 4xl0 20 cm~ 2 , we ob- 
tain a best-fit value for the oxygen abundance of 2.0to'J 
times the solar value (Anders & Grevesse 1989). On the 
other hand, two instrumental absorption features at 23.05 
and 23.35 A contribute ~ 2xl0 17 cm~ 2 to the oxygen col- 
umn density (de Vries et al. 2003) and could also explain 
most of the detected overabundance of oxygen. 

If the emission line at 19.19 A is redshifted O VIII, we 
might expect spectral features from other ionization states 
of oxygen. The signature of He-like oxygen is a triplet of 
X-ray lines at 21.6, 21.8, and 22.1 A. Some residuals are 
in fact present in the RGS spectrum around 21.8 A, which 
corresponds to either the rest-wavelength of the intercom- 
bination line or a resonance line at z=0.01. Fitting the 



4 



Tiengo et al.: Discovery of a redshifted X-ray emission line in the symbiotic neutron star binary 4U 1700+24 



'■■■>v-. >% 



^ o 
3 o - 



CV2 r 
O LTD 



L_l I I I I I I I I I I I I I I I I I I I I L 

10 15 20 25 30 35 

Wavelength (Angstrom) 

Fig. 2. RGS1 (filled squares) and RGS2 (open circles) 
spectra of 4U 1700+24 . The residuals in the middle panel 
are relative to the model derived from the PN data anal- 
ysis (see table 1) multiplied by a factor 0.8 to account for 
the discrepancies in flux estimates between the PN and 
RGS. The lower panel shows the residuals relative to the 
same model plus an emission line at 19.19 A. For clarity, 
the data have been rebinned to a bin size larger than the 
RGS energy resolution. The data gaps at 10.5-14 A for 
RGS1 and 20 - 24 A for RGS2 correspond to the CCDs 
lost due to electronic failure in the early phase of the mis- 
sion. 



residuals at 21.8 A with a Gaussian of the same width 
as the O VIII line candidate (a = 4 eV) gives a 3<r upper 
limit for the equivalent width of 5 eV. Fully ionized oxygen 
should produce an O VIII radiative recombination contin- 
uum at around 14 A (14.2 A for plasma at rest and 14.4 
A at z = 0.01). The RGS spectrum shows some residuals 
at this location. Including a recombination emission edge 
at 14.4 A with a 50 eV width in the fit gives a normaliza- 
tion of (6.1±3.1)xl0 -4 photons cm" 2 s" 1 . The addition 
of this component is statistically significant, although it 
is quite broad and could also result from uncertainties in 
the spectral continuum model. The RGS spectra were also 
rebinned using different criteria, but no other spectral fea- 
tures could be identified. 

Due to the poorer energy resolution of the PN camera, 
no line at 19.19 A is significantly detected in the PN data. 
However, if we fix the energy and the width of the line to 
the values derived from the RGS data, the upper limit to 
the line equivalent width is 12 eV, which is consistent with 
the RGS results. 



4. Discussion 

The high sensitivity of XMM-Newton at low energies and 
the high spectral resolution of the RGS instrument have 
revealed new features in the X-ray spectrum of the un- 
usual interacting binary 4U 1700+24 . Although the PN 
spectrum above 1 keV is consistent with the rather vari- 
able spectra seen in previous observations (Masetti et al. 



2002), we clearly detect both a soft excess and an emission 
line at 19.2 A. 

The best fit to the soft excess was found by adding a 
broad Gaussian line component to the high-energy contin- 
uum model (blackbody + Comptonization) . But the inter- 
pretation of this feature as an emission line is problematic. 
No strong emission lines are expected at the wavelength 
of the soft excess for either z =0.01 or z =0. In addition, 
both the intensity and broadening of this line would be 
much larger than for the line at 19.2 A. 

As mentioned above, we have also discovered an emis- 
sion line at 19.2 A. Since only very low emissivity lines are 
consistent with the observed wavelength of A = 19.19+q 09 
A, the line may be O VIII at redshift ~0.01. For most high- 
emissivity lines close to this wavelength, a much larger 
redshift/blueshift would be required. For example, among 
the other strong lines expected in the spectrum of an ac- 
creting X-ray binary, the O VII and Ne IX triplets are 
the nearest candidates for the line identification, at longer 
and shorter wavelengths, respectively. For O VII, the line 
would be blue-shifted by a factor > 0.1, which, in case of 
Doppler shift, requires that the emitting plasma is mov- 
ing towards us at more than 10% of the speed of light. On 
2002 July 29 radio observations were performed to look for 
a possible jet, but no source was detected with an upper 
limit of 1.0+0.7 mJy at 15 GHz (G. Pooley 2002, private 
communication) . 

On the other hand, the identification of the line at 
19.2 A with the Ne IX triplet would imply redshifts of 
0.40, 0.41, and 0.42 for the forbidden, intercombination 
and resonance lines, respectively. These values could be 
interpreted as gravitational redshift from close to the sur- 
face of a neutron star, as they are consistent with most 
of the equations of state for neutron stars composed of 
normal nuclear matter, and they are just slightly larger 
than the redshift z=0.35 found by Cottam et. al. (2002) in 
the spectral analysis of X-ray bursts from EXO 0748-676. 
However, due to the lack of identifications of other spec- 
tral features at z ~ 0.4 or indication for overabundance of 
neon, we consider O VIII to be a better candidate for the 
observed emission line. 

The relatively small redshift of the O VIII line (z ~ 
0.01) can be interpreted in several ways, and we briefly dis- 
cuss two possible scenarios. A value of z = 0.008-0.014 cor- 
responds to the gravitational redshift of a photon emitted 
at a distance of 35-60 Schwarzschild radii from a compact 
object: this interpretation would exclude the possibility 
that 4U 1700+24 is a white dwarf and would correspond 
to a distance of 1.5-2.6xl0 7 cm from a 1.4 M neu- 
tron star. In a photoionized plasma, assuming that most 
of the oxygen is H-like, the measured O VIII line lumi- 
nosity of ~10 40 photons s _1 gives an emission measure of 
EM ~ 6x 10 53 cm" 3 (taking an abundance of ~1.6xl0" 3 , 
as derived from the oxygen edge fit and assuming that the 
overabundance is intrinsic to the X-ray source) . Assuming 
a spherically symmetric geometry, if the O VIII line is 
emitted at ~ 2xl0 7 cm from the neutron star, we can es- 
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timate a density of n ~ 4 x 10 15 cm~ 3 . For the measured 
continuum luminosity of 10 34 erg s _1 , the ionization pa- 
rameter is £ =L cont /(n R 2 ) ~5000 erg cm s" 1 . For such 
a high ionization parameter, most of the oxygen should 
be fully ionized and therefore we would not expect to see 
a prominent O VIII emission line. O VIII dominates the 
ionization states of oxygen for £ ~ 100 erg cm s _1 , which 
means that either we have used an oversimplified model 
(assuming, for example, that the line emission region is 
symmetric and homogeneous) or that the O VIII line is 
emitted at larger distance from the central X-ray source 
and therefore that the redshift is not solely gravitational. 

In an alternative scenario, the same redshift can be 
produced by Doppler effects if the emitting plasma is mov- 
ing away from us at a speed of 2000-4000 km/s. This 
velocity is about two orders of magnitude larger than 
the wind velocity of an M-type giant (Reimers 1977), 
as well as the proper motion and possible orbital veloc- 
ity of the binary system (Galloway, Sokoloski & Kenyon 
2002). Although the high luminosity, hard X-ray spec- 
trum, and rapid X-ray variability of 4U 1700+24 make 
it rather unlikely, the present data do not rule out the 
possibility that the accreting object is a white dwarf. The 
detection of Dopplcr-shiftcd lines of highly ionized ele- 
ments has been reported for some supersoft sources and 
interpreted as the signature of collimated outflows ("jets") 
coming from the accreting white dwarf (e.g., Cowley et al. 
1998). Furthermore, the jet velocity is typically similar 
to the value we derive for the redshifted O VIII line. In 
our case, however, no corresponding blueshifted line is de- 
tected, which implies either a unipolar jet or some special 
geometry to obscure the approaching jet. 

For a neutron star with a 10 12 gauss magnetic field and 
a luminosity of 10 erg s _1 , the magnetospheric radius is 
~ 3x 10 9 cm (Hayakawa 1985). Although X-ray pulsations 
have not been found from 4U 1700+24 , the fact that it 
is probably viewed close to face-on (Galloway et al. 2002) 
means that pulsations would be undetectable even if it is 
an accreting neutron star with high magnetic field as the 
line of sight is almost aligned with the neutron star spin 
axis. At a distance of ~ 3x 10 9 cm from a 1.4 M Q neu- 
tron star, a free particle would move radially towards the 
neutron star with a velocity of ~ 3400 km/s. Therefore, 
the redshift we measure could originate from close to the 
magnetospheric radius, where the plasma density is in- 
creased by the funneling effect produced by the magnetic 
field lines. Repeating the estimates of the ionization pa- 
rameter described above, at a distance of 3x 10 9 cm the 
corresponding density is n ~ 2 x 10 12 cm~ 3 and £ ~ 500. 
These two values are fairly consistent with the presence of 
an O VIII line. 

Redshifted emission lines have also been reported in 
Cygnus X-3. Paerels et al. (2000) reported a redshift cor- 
responding to 750-800 km/s for all the emission lines de- 
tected in a Chandra HETGS observation of Cygnus X- 
3. In two more recent Chandra observations, a slightly 
smaller redshift of 270-550 km/s, in addition to a Doppler 
modulation of about 150 km/s, was found (Stark & Saia 



2003). The second effect was interpreted as being due to 
the binary orbital motion, but no convincing explanation 
was given for the global redshift. Although the X-ray lumi- 
nosities and observed redshifts of 4U 1700+24 and Cygnus 
X-3 are rather different, a common mechanism could pro- 
duce the line redshifts. For both 4U 1700+24 and Cygnus 
X-3, a low orbital inclination has been proposed (Galloway 
et al. 2002; Mioduszewski et al. 2001). Thus, Doppler 
shifts due to the bulk velocity of matter receding from 
us as it falls onto the compact object could be the origin 
of the redshifted lines. 
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